# reporting unit locationsdat<-read_csv("/Users/jeffbaldock/Library/CloudStorage/GoogleDrive-jbaldock@uwyo.edu/Shared drives/wyo-coop-baldock/UWyoming/Snake River Cutthroat/Methods/GSI Sampling/Snake_GSI_field data_2020-2022_GenOnly_DropNoDrop_090823edit.csv")%>%filter(collection.type=="baseline")%>%group_by(site)%>%summarize(lat =unique(start.lat), lon =unique(start.lon))%>%ungroup()%>%mutate(repunit =site)# define reporting groups according to GSIdat$repunit[which(dat$repunit=="bluecrane_ford")]<-"cody_bluecrane"dat$repunit[which(dat$repunit=="cody_moment")]<-"cody_bluecrane"dat$repunit[which(dat$repunit=="cowboycabin_springchannels")]<-"cowboycabin_NA"dat$repunit[which(dat$repunit=="threechannel_eastfork")]<-"threechannel_NA"dat$repunit[which(dat$repunit=="threechannel_westfork")]<-"threechannel_NA"dat$repunit[which(dat$repunit=="fish_lower")]<-"fish_NA"dat$repunit[which(dat$repunit=="fish_upper")]<-"fish_NA"dat$repunit[which(dat$repunit=="slate_lower")]<-"slate_NA"dat$repunit[which(dat$repunit=="slate_upper")]<-"slate_NA"dat$repunit[which(dat$repunit=="spread_northfork")]<-"spreadnf_flagstaff"dat$repunit[which(dat$repunit=="flagstaff_NA")]<-"spreadnf_flagstaff"dat$repunit[which(dat$repunit=="fall_upper")]<-"fall_coburn"dat$repunit[which(dat$repunit=="fall_lower")]<-"fall_coburn"dat$repunit[which(dat$repunit=="coburn_NA")]<-"fall_coburn"dat$repunit[which(dat$repunit=="littlegreys_lower")]<-"littlegreys_steer"dat$repunit[which(dat$repunit=="littlegreys_upper")]<-"littlegreys_steer"dat$repunit[which(dat$repunit=="steer_NA")]<-"littlegreys_steer"dat$repunit[which(dat$repunit=="pacific_lower")]<-"pacific_NA"dat$repunit[which(dat$repunit=="pacific_upper")]<-"pacific_NA"dat$repunit[which(dat$repunit=="deadman _greys")]<-"deadman_greys"# drop reporting groups that cannot be easily combineddat<-dat%>%filter(!repunit%in%c("schwabacher_NA", "southbuffalofork_NA", "jack_NA", "bacon_NA", "bear_NA", "sheep_NA"))# filter to most downstream sitesdat<-dat%>%mutate(ds =ifelse(row_number()%in%c(3,4,7,13,18,24,25,29,31,33,34,16,38,37,17,40,44,41,42,80,45,48,51,53,54,58,61,63,64,66,90,71,72,74,77,78,82,84,86,87,91,94), 0, 1))dat<-dat%>%filter(ds==1)%>%select(repunit, lat, lon)%>%group_by(repunit)%>%summarize(lat =mean(lat), lon =mean(lon))dat<-dat[-33,]# drop lake crk# edit lat/longdat$lat[which(dat$repunit=="cottonwood_grosventre")]<-43.55226dat$lon[which(dat$repunit=="cottonwood_grosventre")]<--110.25895dat$lat[which(dat$repunit=="threechannel_NA")]<-43.55138dat$lon[which(dat$repunit=="threechannel_NA")]<--110.79136dat$lat[which(dat$repunit=="cody_bluecrane")]<-43.43554dat$lon[which(dat$repunit=="cody_bluecrane")]<--110.82925dat$lat[which(dat$repunit=="flagstaff_NA")]<-43.78173dat$lon[which(dat$repunit=="flagstaff_NA")]<--110.28081dat$lat[which(dat$repunit=="dell_NA")]<-43.23103dat$lon[which(dat$repunit=="dell_NA")]<--110.42322dat$lat[which(dat$repunit=="lowerbarbc_NA")]<-43.54686dat$lon[which(dat$repunit=="lowerbarbc_NA")]<--110.78555dat$lat[which(dat$repunit=="spread_southfork")]<-43.74837dat$lon[which(dat$repunit=="spread_southfork")]<--110.31922dat$lat[which(dat$repunit=="blackrock_lower")]<-43.82523dat$lon[which(dat$repunit=="blackrock_lower")]<--110.35208# create spatial data object and project to crsdatsp<-vect(dat, c("lon", "lat"), crs ="+proj=longlat")datsp<-terra::project(datsp, basin)# write out reporting unit locationswrite_csv(dat, "/Users/jeffbaldock/Library/CloudStorage/GoogleDrive-jbaldock@uwyo.edu/Shared drives/wyo-coop-baldock/UWyoming/Snake River Cutthroat/Analyses/Snake River GSI Quarto/Landscape Covariates/Location Data/RepUnit_LatLong.csv")writeVector(datsp, "/Users/jeffbaldock/Library/CloudStorage/GoogleDrive-jbaldock@uwyo.edu/Shared drives/wyo-coop-baldock/UWyoming/Snake River Cutthroat/Analyses/Snake River GSI Quarto/Landscape Covariates/Location Data/RepUnit_SpatialLocations.shp", overwrite =TRUE)
Code
datsp<-vect("/Users/jeffbaldock/Library/CloudStorage/GoogleDrive-jbaldock@uwyo.edu/Shared drives/wyo-coop-baldock/UWyoming/Snake River Cutthroat/Analyses/Snake River GSI Quarto/Landscape Covariates/Location Data/RepUnit_SpatialLocations.shp")